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Abstract: The modeling ability of a stand-alone version of the Simple 
Biosphere Model 2 (SiB2) was tested mainly through diagnosing the 
simulated latent heat (LE), sensible heat (H), CO 2 flux, and air tempera¬ 
ture at the Tongyu field observation station (44°25N, 122°52'E, 184 m 
elevation) of Coordinated Enhanced Observing Period (CEOP), where 
the land cover is cropland and grassland. In the whole year of 2003, the 
canopy height and the leaf area index was variable. During non-growth 
period, the surface would become bare, while during the growth period, 
the canopy height could reach 2.0 m high over cropland and 0.8 m high 
over grassland, respectively, and max leaf area index could reach 4.2 and 
2.4, respectively. The model was initialized with measurement and driven 
by half-hourly atmospheric observations. The simulation values for 2003 
were compared against measurements. Results show that the model is of 
a good ability of simulating the hourly latent heat (LE), sensible heat (H), 
C0 2 flux and temperature during the growth period. Moreover, the daily 
LE, H and C0 2 flux simulated by SiB2 could reflect their yearly change 
reasonably. However, the model may overestimate the H generally. 
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Introduction 

Grasslands occupy a significant portion of the earth's land area. 
Grassland ecosystems are diverse and have evolved in response 
to large seasonal and annual fluctuations in moisture, tempera¬ 
ture, and grazing intensity (Knapp et al. 1998; Schuman et al. 
2000). In recent decades, severe land degradation and desertifi¬ 
cation have occurred in China due to overgrazing, excessive 
mowing, and particularly crop cultivation, which are caused by 
the food demands for an increasing human population (Zhou 
2002). The drought cropland and grassland in monsoonal Asia 
play an important role in the East Asia Monsoon energy balance 
(Fu et al 1999; Fu 2002; 2003). Degraded land provide a greater 
potential for sequester carbon. The improvement of degraded 
grasslands provides economic and sociological benefits as well. 
Tongyu flux station, a station in Coordinated Enhanced Observ¬ 
ing Period (CEOP) (www.ceop.net) and founded in October 
2002, was established for a long-term, continuous observation of 
the human activities and the exchange among surface atmosphere, 
lower-level atmosphere and biogeocenose. The constant data are 
very useful for understanding the impact of land use and water 
use on regional energy circulation and water circulation. The 
data obtained from this station have been used successfully to 
test the simulation abilities of several well known land surface 
models and ecosystem models, including IBIS model (Jiang et al. 
2008), and BIOME-BGC model (Wang et al. 2006). 

Atmospheric general circulation models are used for climate 
simulation and weather forecasting, requiring the fluxes of radia¬ 
tion, heat, water vapor, and momentum across the 
land-atmosphere interface that has been specified (Sellers et al. 
1997). Land use and land cover (LULC) is the most important 
factor affecting the regional and global climate (Pielke et al. 
1990; Pitman et al. 1999). The interaction between the vegetation 
cover and the lower-level atmosphere plays an important role in 
regulating the atmospheric circulation and then the various cli¬ 
mate factors (Dickson et al. 1993; Irannejad et al. 1998). 

The atmosphere, vegetation, and soil system are dynamically 
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coupled through the physical processes, which produce the 
transport of thermal energy and mass across the interface. 
Therefore, modeling based on physic is an important tool for 
studying the coupled system (Sridhar et al. 2002). In this study, 
the micrometeorological data at Tongyu flux station located in 
Northeast China was used to simulate latent heat, sensible heat 
and C0 2 with a stand-alone version 2 of the Simple Biosphere 
Model (SiB2) (Sellers et al. 1996a; 1996b). The SiB2 was origi¬ 
nally designed to work with atmospheric general circulation 
models (GCMs). Randall et al. (1996) compared the results of a 
simulation of SiB2, which had been coupled to a GCM, against 
the results of a controlled simulation of SiB2, which had been 
coupled to a bucket surface hydrology model. They found that 
SiB2 produced warmer and drier surface and atmospheric 
boundary layer in the former work than the controlled work, 
resulting in increased surface sensible heat and decreased latent 
heat over the continents. Gao et al. (2004) used an off-line SiB2 
to simulate daily and seasonal variations of the components of 
the surface energy balance and the C0 2 flux, though underesti¬ 
mating the sensible heat flux and latent heat flux by 10% and 
11%, respectively. Gao et al. (2004) found that SiB2 could esti¬ 
mate soil wetness reasonably in the Tibetan prairie and still gen¬ 
erated warmer (colder) ground surface in day-time (night-time) 
when they underestimate net radiation by 11%, and overesti¬ 
mated sensible, latent, and soil heat fluxes by 8%, 3%, and 11%, 
respectively. 

In order to improve the current understanding of heat fluxes, 
and C0 2 exchange in semi-arid cropland and degraded grassland, 
we had designed an intensive cropland and grassland experiment 
in China’s Jilin province since October 2003. This site is a typi¬ 
cal semi-arid land in east-north china. Two different kinds of 
land covers, cropland and degraded grassland were considered to 
understand the characteristic of the semi-arid land. 

Materials and methods 

Study area 

Measurements at Tongyu field observation station (44°25'N, 
122°52'E, 184 m elevation), one of the referenced sites of CEOP, 
had been made since October 2002. This field occupied cropland 
with maize planted and degraded grassland either, and the sites in 
this field were flat, with approximate 5 km distance between any 
two sites. Soil at the experimented site was sandy loam, contain¬ 
ing 60% sand, 35% clay, and 5% silt, respectively. Over the 
cropland, maize was cultivated from May to earlier October 
every year, so during the growth period the canopy height could 
be increased from 0.0 to 2.0 m, but during the other times the 
land cover was bare soil surface; over the degraded grassland, 
the grass maximum height could only reach 0.08 m. 

Micrometeorological measurements 

At every site, a 20-m meteorological observation tower was set 
up, and on the tower, wind speed (by Met one, 034 A-C), tem- 
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perature and humidity (by Vaisal, Hmp45C-L) were measured 
over five layers, and the wind direction (17 m, by Metone, 
014A-L) was measured over 17 m layer. On the ground, the sur¬ 
face temperature (by temperature Hitester, APOGEE, IRTSO-P), 
air pressure and precipitation were measured. Upward and 
downward solar radiations (by CM21, Kipp & ZONEN) and 
upward and downward long wave radiations (by CG4, Kipp & 
ZONEN) were measured at 2.0 m and 3.0 m, respectively. Six 
layers soil temperature (2, 5, 10, 20, 50, 80 cm, by STPO1-L50, 
107-L) and five layers soil volumetric Water Content (5, 10, 20, 
40, 80 cm by CS616-L) were measured. Soil heat fluxes (by 
HFP01SC-L50) were measured at the soil depths of 5 and 10 cm. 
The sample frequency for those variations was 0.5 Hz. The eddy 
covariance instruments including a three-dimensional sonic 
anemometer (CSAT3), and humidity and C0 2 pulsation instru¬ 
ment (LI-COR, CS7500) were set at 3.0 m over cropland and at 
2.0 m at grassland, respectively. The sample frequency was 10 
Hz. 

The micrometeorological measurements and fluxes were in¬ 
vestigated during CEOP-EOP3 period (mainly in 2003). 

SiB2 model 

Overview 

The SiB2 is a simple but realistic biosphere model for calculating 
the transfer of energy, mass and momentum between the atmos¬ 
phere and the vegetated surface of the earth. It is designed for 
use in atmospheric general circulation models (Sellers et al. 
1986). The detailed structure of SiB2 is described by Sellers et al. 
(1996a). The SiB2 uses two-stream approximation for modeling 
radiation. The SiB2 incorporate a realistic canopy photosynthe¬ 
sis-conductance model to describe the simultaneous transfer of 
C0 2 and water vapor into and out of the vegetation, respectively. 

Forcing data 

In the off-line Simple Biosphere Model 2, 10 meteorological 
forcing variables, including downward short-wave radiation, 
downward long-wave radiation, net radiation, vapor pressure, air 
temperature, wind speed, precipitation, air pressure, soil surface 
temperature, and soil deep temperature, were required. These 
forcing data, except for the precipitation, could be got basically 
by half-hourly measurements at Tongyu flux station. The pre¬ 
cipitation could be provided on a day basis after July 2003. To 
meet the demand of the model, we incorporate the half-hourly 
vapor pressure and the daily precipitation to get the half-hourly 
precipitation as shown in Fig. 1, through assuming that the 
half-hourly precipitation was proportionate to the half-hourly 
vapor pressure daily because of lack of hourly measurement. 
Three layers (5 cm, 10 cm, and 20 cm) of half-hourly soil volu¬ 
metric water content were also required. The reference height 
was determined to be 2.95 m high over cropland and 1.35 m high 
over grassland. 

Parameter setting 

SiB2 required users to input the soil and vegetation types. The 
soil at cropland and grassland contained mainly sand and loams, 
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therefore we selected the “sandy loam” type soil from Sellers et 
al. (1996b, in the Table 4). Over the cropland, maize was culti¬ 
vated, but the degraded grassland was dominated by C 3 
short-grass species. Therefore, together, the category was defined 
as “agriculture or C 3 grassland” (biome type 9) (Sellers et al. 
1996b, in the Table 2). We obtained the time-invariant vegetation 
parameters from Sellers et al. (1996b, in the Table 5). The bulk 
boundary layer resistance coefficient Q and the ground to can¬ 
opy air-space resistance coefficient C 2 were set to 15.1 and 131.0 
for cropland, and 4.0 and 70.0 for grassland, respectively (Sellers 
et al. 1996b). During the whole year, the canopy heights changed 
from the minimum to the maximum, and then to the minimum 
again. The leaf area index, canopy greenness fraction, canopy 
cover fraction, canopy roughness length and canopy zero plane 
displacement also had changes. Considering the monthly mean 
leaf-area index data from the MODIS 
(ftp://primavera.bu.edu/pub/datasets), we combined the discon¬ 
tinuous observation with the logistic model (Lin et al. 2003), and 
then gave the variation of LAI in 2003 as shown in Fig. 2. The 
LAI had an abrupt decrease due to the harvest of maize. The 
canopy heights over the cropland and degraded grassland were 
shown in Fig. 3. We set the roughness length Z 0 =Z 2 /7.6 (Tanner 
et al. 1960), the canopy zero plane displacement D=0.7Z 2 (Stan- 
hill 1969; Thom 1971), the canopy cover fraction (V) to 0.8 over 
cropland and 0.6 over grassland. The canopy greenness fraction 
(N) equaled to 1 during the growth period and 0 during the other 
times. 
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Fig. 2 Yearly mean leaf-area index (LAI) data over cropland and 
degraded grassland. 


grassland. Because the flux changed mainly in growth period, we 
selected 213-233 Julia days of the year, in which period, maize 
grew to mature over the cropland and its leaf area index could 
reach 4.2. When maize matured, its respiration decreased. Fig. 4 
showed the time series of the latent heat (LE), sensible heat (H), 
C0 2 flux and air temperature (T) modeled by SiB2 or obtained 
by the direct measurements. It was obvious that SiB2 had gener¬ 
ated consistent temporal variations of LE, H, C0 2 flux and T. 
The squared correlation coefficients were shown in Table 1. The 
scatterplots of the modeled LE, H, C0 2 flux and T, and the direct 
measurements were given in Fig. 5. It was apparent that SiB2 
overestimated LE in nighttime and H in all day. Moreover, it also 
estimated more C0 2 flux released by maize and soil land cover 
in nighttime. Nevertheless, it estimated the air temperature per¬ 
fectly. 
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Fig. 3 Heights of canopy top (Z 2 ) and canopy bottom (Zi) over crop¬ 
land and degraded grassland, respectively. 
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Fig. 4 Hourly change (213-233 th days) of latent heat flux (LE), sensi¬ 
ble heat flux (H), and C0 2 flux modeled using SiB2 against direct 
measurements over cropland. 


Results 

Hourly modeling of the growth period 

Based on the half-hourly measurements, the half-hourly steps 
were used to simulate the annual changes over cropland and 


Table 1. Squared correlation coefficients of daily latent heat (LE), 
sensible heat (H), C0 2 flux and temperature (T) in cropland and 
grassland for growth season. 


Land use types 

LE 

H 

C0 2 flux 

T 

cropland 

0.8977 

0.809 

0.864 

0.9991 

grassland 

0.8911 

0.8664 

0.8063 

0.9971 
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Fig. 5 Scatterplots of latent heat flux (LE), sensible heat flux (H), and 
CO 2 flux modeled using SiB2 against direct measurements over crop¬ 
land. 
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Fig. 6 Hourly change (213-233th days) of latent heat flux (LE), sensi¬ 
ble heat flux (H), and C0 2 flux modeled using SiB2 against direct 
measurements over grassland. 


Fig. 6 showed the comparison between the modeled and the 
direct measurements over grassland, and Fig 7 was its scatter¬ 
plots. SiB2 overestimated the LE and the H. Although it overes¬ 
timated the air temperature, its squared correlation coefficient 
reached 0.9971 as shown in Table 1. 

Hourly modeling of mean day 

The ability of hourly simulation in a day was investigated by 
comparing the modeled half hourly latent heat (LE), sensible 
heat (H), C0 2 flux and air temperature (T) with the several-day 
averaged value of the direct measurements. Figs. 8-9 showed that 
SiB2 had a good hourly simulating ability for cropland and 
grassland, and their squared correlation coefficients were shown 
in Table 2. 


Table 2. Squared correlation coefficients of hourly latent heat (LE), 
sensible heat (H), C0 2 flux and temperature in cropland and grass¬ 
land for growth season. 


Land use types 

LE 

H 

C0 2 flux 

T 

cropland 

0.9857 

0.9521 

0.9634 

0.9997 

grassland 

0.9907 

0.9879 

0.9703 

0.9997 
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Fig. 7 Scatterplots of latent heat flux (LE), sensible heat flux (H), and 
C0 2 flux modeled using SiB2 against direct measurements over 
grassland. 
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Fig. 8 Hourly change of a mean day over cropland (averaged from 
DOY (Julia days of the year) 200 to DOY 210). 
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Fig. 8 showed the mean variations of DOY (Julia days of the 
year) 200-210, in which period, the maize grew rapidly with the 
strongest respiration, and its leaf area index could reach 4.2. This 
proves that SiB2 had done a good hourly simulation. It also had 
high squared correlation coefficients over cropland except for 
overestimating H (Table 2). In the nighttime, there was no pho¬ 
tosynthesis of maize and no heat by radiation either, so the LE, 
the H, and C0 2 absorbed by maize and soil decreased to the 
lowest value in a day. When sun rose, the temperature near the 
surface rose and the photosynthesis began, so the LE and the H 
increased (reach their maximum at noon), and so did the C0 2 
absorbed by the surface. However, at noon, the maize C0 2 ab¬ 
sorption decreased sharply due to high temperature caused 
stomatal closure. All the variables would descend with the de¬ 
clining radiation in the afternoon. SiB2 could simulate all these 
process reasonably. SiB2 could also simulate this phenomenon 
over grassland. 

Daily modeling 

The half-hourly time steps were used to investigate the ability of 
daily simulation of the model. During a year, the characteristics 
of the land surface often changed with the variety of the radiation 
and the anthropic activities. Over the cropland, the maize was 
planted in May and harvested in early October, and during the 
other times the land was bare soil. The measurements in the first 
45 days from 10 April to 15 May were discarded for its lack of 
measure and strongly anthropic activities. It is obvious that SiB2 
consistently generated temporal variations of the LE, the H, the 
C0 2 , and the T (Fig. 10). The squared correlation coefficients are 
shown in Table 3. The same variations were also shown on 
grassland (Fig. 11). SiB2 gave a perfect simulation of the air 
temperature. During late May and June, the rapid growth of the 
maize increased the model’s disturbance. SiB2 also gave a rea¬ 
sonable estimate of the C0 2 flux, and the minus values represent 
the C0 2 absorbed by plants and soils. However, the measurement 
error in the non-growth period might occur due to less C0 2 ex¬ 
change caused by less vegetation cover and low temperature. 


Table 3. Squared correlation coefficients of daily latent heat (LE), 
sensible heat (H), C0 2 flux, and temperature in cropland and grass¬ 
land for whole year. 


Land use types 

LE 

H 

C0 2 flux 

T 

cropland 

0.8535 

0.7353 

0.7542 

0.9998 

grassland 

0.9184 

0.8539 

0.5903 

0.9981 
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Fig. 11 Comparison of latent heat flux (LE), sensible heat flux (H), 
C0 2 flux, and air temperature modeled using SiB2 against direct 
measurements over grassland. 

The simulated and observed values had some obvious differ¬ 
ence at some moment for the four parameters. For example, Fig. 
10 showed a poor LE simulation. This might be caused by com¬ 
plex reasons such as the observation errors (almost all flux tower 
measurement in fact exists a problem: energy budget is not close, 
thus LE and H may have a gap to actual flux), the model limita¬ 
tion (some physical processes may be simplified incorrectly), 
and human activity effects (including possible irritation, fertili¬ 
zation and pasture, etc.) 
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Fig. 10 Comparison of mean daily (24-hour) latent heat flux (LE), 
sensible heat flux (H), C0 2 flux, and air temperature (T) modeled 
using SiB2 against direct measurements over cropland. 


Summary and conclusions 

In this study, the modeled results of the energy flux and C0 2 flux 
were obtained over a cropland and degraded grassland in 
north-east China in 2003. The daily patterns of the SiB2 modeled 
latent heat flux (LE), the sensible heat flux (H), the C0 2 flux and 
the air temperature were consistent with the direct measurements 
during the growth period, especially with the day-averaged value. 
SiB2 showed a good simulation ability of the hourly changes, 
and also gave a reasonable simulation of the seasonal varieties in 
a whole year. 

However, many gaps between modeling and observations in¬ 
evitably exist. Besides observation methods and condition are 
imperfect, the model simplification is another key factor. For 
example, LE simulation can be closely correlated with soil heat 
process, which was approximately computed using forc¬ 
ing-restore method. In fact, this only is a first order approximate 
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and not suitable for real non-linear process. Hence, we think that 
soil thermo process and vegetation growth and physiological 
process should be improved for more realistic simulation of the 
water flux with heat flux in the future. 
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